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Abstract 

The observed J/tfj suppression in Pb + Pb collisions shows a drop at those large values of 
the transverse energy Et which arise from fluctuations. The validity of existing models for J/ip 
suppression can be extended into this domain of Et by introducing an ad hoc factor proportional 
to Et- We propose a formalism in which the influence of Et fluctuations on the J/ip suppression 
can be calculated and discuss the conditions under which the ad hoc factor is obtained. 
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The observed dependence on transverse energy of J/ip production in Pb+Pb collisions at 158 A 

O ■ GeV is the first unambiguous signal for an anomalous mechanism for charmonium suppression, i.e. 
O ' 

one which goes beyond what is already observed in proton-nucleus collisions and in reactions with 



light ions. In order to identify the detailed nature of the anomalous mechanism, it is important to 



i 

O : investigate each of its features using a mi— of adjustable parameters. The past experience has 
Ch shown that models succeed to reproduce the data only after several parameters are adjusted. Three 
' years ago, the NA50 collaboration has re-measured the J/ip suppression in Pb+Pb collisions in the 
'X ; region of large transverse energy E T and has corrected its earlier result Q. Rather than being flat 
5r as a function of Et the new data show a drop at Et ~ 100 GeV. For an overview of the field, the 
data and their interpretation we refer to the reviews [0, || and the proceedings of the Quark Matter 
'99 conference [§]]. 

Within the picture of the quark-gluon plasma the newly observed drop is interpreted as the 
expected onset of J / ip melting . Two other groups, Capella et al. and Blaizot et al. Q7p , modify 
their previously successfully proposed expressions for the observed J/ip suppression by introducing 
the factor 

at the appropriate place, without the need of a new parameter. They argue that the newly discovered 

drop for J /if; suppression arises in the regime where the increase in Et is not caused by a decrease in 

1 Present address: Physik Department, TU Miinchen, James-Franck-StraBe, D-85747 Garching, Germany, e-mail: 
polleri@physik.tu-muenchen.de 



1 



impact parameter b, but rather by fluctuations in Et around the mean value -&r(&), which is deter- 
mined by the collision geometry. Larger values of Et lead to a correspondingly larger suppression. 
For the comover description by Capella et al. the modification by the factor e has a small effect 
and fails to describe the dataQ, while the same modification to the cut-off model by Blaizot et al. 
quantitatively describes the newly observed break. These partial and full successes, respectively, call 
for a derivation of the ad hoc formula. This is what is attempted in this paper. 

The expression for the charmonium (ip) production cross section in A + B collisions is usually 
written as 

da\ B 



dEn 



where 



-(E T ) = ^ p Jd 2 bP T (E T ,b) (2) 
x J d 2 s T A T B S nud (b, s) S FS i(b, s) , 

1-exp (-a$T A ) l-exp(-gffir B ) 

bnud abs rp a bs rp (.-J J 

represents the part of the suppression (also present in pA collisions) which is related to the propa- 
gation of the charmonium through both nuclei, the ip being destroyed in ipN collisions with a rate 
determined by the absorption cross section (X^y-- ^Ve denote by T A — T A (s) and T B — T B (b — s) the 
nuclear thickness functions. Pt represents the probability of observation of Et in events with impact 
parameter b, and is normalized to 1 when integrated over Et- The anomalous part Sfsi accounts 
for final state interactions with the produced quarks and gluons (QGP) or the hadrons (comovers). 
In the comover approach |J one writes 

S c F ° SI (b, s) = exp{-a co N c y °(b, s) \n(N c y °(b, s)/N f )}. (4) 

The suppression function depends on the density of comovers N£° at the rapidity y of the ip, while a co 
describes the ip absorption by a comoving meson and is usually taken as the adjustable parameter of 
the theory (of the order of 1 mb). The corresponding density of comovers in pp collisions is denoted 
by Nf. 

In the cut-off model one assumes 

S^j(b,s) = Q(n;-n p (b,s)), (5) 
where the (mean) density of participant nucleons in impact parameter space is 

n p (b, s) = T A [1 - exp(-a pp T B )} + (T A <-> T B ) . (6) 

2 In a recent paper Capella et al. || argue, that the drop shown in the data is misleading, since the data are 
evaluated with respect to minimum bias events, whose Et distribution differs from the Et distribution of events in 
which also a ip is observed. 
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The "critical participant density" n p , a parameter adjusted to n\ 



p 



c 



3.7 fm" 2 , 



represents the density 



above which if) absorption is 100 % effective. The QGP phase transition may (but must not) be the 
mechanism for the critical transition. In ref. @ the B-function in eq. @ is smeared at the expense 
of a further parameter A, which is obtained from a fit to the data: 



In order to describe the endpoint behavior of charmonium suppression at large Et, the factor e 
of eq. (P is introduced in front of the factor Ny° in eq. @), with the argument that the number of 
comovers fluctuates proportionally to the observed Et, a very plausible assumption. In eq. fl5|) for 
the cut-off model, the factor e multiplies n p , the mean number of participants. This modification is 
not obvious. Why should the number of participants in the tube, where a if) is produced, fluctuate 
proportionally to the global value of E? ? For instance, in a central collision practically all nucleons 
participate. Their number equals A + B and cannot fluctuate. However, the number of produced 
particles or the produced energy density does fluctuate. Therefore, if the introduction of the factor 
e in front of n p has to make sense, we must interpret n p and n c p as being proportional to the energy 
or particle density in the tube where J /if) is produced. This is indeed the point of view of ref. 0. A 
proportionality factor is irrelevant in the 0-function, but has bearing on the parameter A. 

The actual situation is somewhat more complicated: the hadrons which are observed as transverse 
energy are measured in a pseudo-rapidity interval 1 < rj < 2.3, while the if) is measured in the rapidity 
interval 3 < y < 4. (Recall that at the SPS mid-rapidity corresponds to y ~ 3.) It is not obvious 
to which degree fluctuations in Et in one rapidity interval should have a bearing on the comovers 
or the energy density which suppress the charmonium in another rapidity interval. There must be a 
mechanism for "cross talk" . 

The paper addresses the following two issues:^ 

(i) To establish a framework in which the influence of fluctuations in Et on if> suppression can 
be calculated, to derive a correction factor containing Et and to investigate the limits under 
which it reduces to the ad hoc expression of eq. ([!]). 

(ii) To investigate the importance of "cross talk" , i.e. to calculate the correlation function for the 
coincidence of hadrons in different rapidity intervals. 

To deal with issue (i) we consider a nucleus-nucleus collision at fixed impact parameter b. The 

individual NN collisions produce many particles, possibly first as strings or partons which then 

convert into the observed hadrons. In what follows we will speak of particles, leaving open whether 

3 After the first version of our paper had appeared on the web, Chaudhuri posted a paper addressing similar 
issues. 
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we deal with partons or hadrons. We consider the 6-dimensional phase-space for the particles and 
divide it into cells numbered as i = 0, k, where the index is reserved for the cell in which ip 
is found together with its comoving particles. We denote by rrii the number of particles in the i-th 
cell and by Pi(rrii) their probability distribution. It will be characterized by a mean value and a 
variance of. The probability to find mo particles in cell in an event with a total of M produced 
particles is then 

k 

(p(M,m ) = Y^ Po(m ) . . .p k (rrik)6(M -^2m k ). (8) 

mi,...,m k i=0 

Note that we do not sum over m . We also introduce the probability distribution 

MEt,M) = \ exp[-(E T -e T M) 2 /2d 2 ET ], (9) 

which describes the correlation between the observed transverse energy Et and the total number M 
of particles, ct being the mean transverse energy which each particle contributes to the total Et- 

The correlation between the observed transverse energy Et and the number m of particles, which 
suppress ip is then given by 

P(E T , m ) = £ MEt, M) <p(M, mo) ■ (10) 

M 

If Gaussian distributions are used for the correlation functions <pT and <p and if the sums over M and 
m are replaced by integrals, eq. fllPl) can be evaluated exactly, but a rather cumbersome expression 
results. For the case of many cells, k ^> 1, the result, can be factorized as 

P(£ T ,m ) = P T (E T ,b) P c (m ,m ,£ T ). (11) 

Here the distribution of the fluctuations in Et is 

Pt{EtM = — L= exp[-(E T - E T (b)) 2 /2a 2 ET l (12) 
V 27r(J iT 

where Et is the mean value of the produced transverse energy for a given value of the impact 



parameter. In eq. (12) Pt does not depend on the number m of particles in cell because they 
contribute a negligible amount. However, the distribution function P c of particle number m does 
depend on Et via 

P s {m ,mo,E T ) = f= (13) 

x exp[ -(mo - m (b, s)e(E T )) 2 /2a 2 ] , 
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where the dependence on transverse energy is contained in the expression 



, (Br) . 1+ao (^|M, (14) 

which is close in shape to the one of eq. ([]]). The factor 

0"o e T E T 

a = — (15) 

mo °e t 

depends on the mean value m and the variance o§ of the probability distributions po for the number 
of particles in cell and on the corresponding quantities E (b) and cr ET (b) for the distribution of the 
observed transverse energy. When both distributions are normal, i.e. m = ctq, <j\ = ctEt, each 
of the two factors in eq. ([15]) equals 1 . Then a = 1 and 

e(E T ) = e(E T ) = Sf (16) 

as assumed in refs. || [7[]. A value for the second factor eTE T /a Er in a can be deduced from a fit of 
calculated Et distributions (using eq. (|T2"D) to the measured one. One finds values between 1 and 0.7 
H [Tj. Although we have no direct information on the ratio a^/mo, the multiplicity distribution of 
produced particles in pp collisions are known to be negative binomials Jl0| for which cx^/mo = 1 + n?, 
where k, the order of the binomial, is found to be of order 3 or 4 for pp collisions at ISR energies. 
Therefore «o is a product of two factors, one larger than 1, the other smaller than 1. Although they 
may compensate each other to a large extent, we are not in a position to give a reliable value for ao- 
The derivation of the expression €(Et), eq. (|14"D is the first result of our paper. 

Now we address issue (ii), the effect of incomplete cross talk. In the above derivation it has been 
tacitly assumed that produced particles in cell 0, which contribute to if) suppression, also contribute 
to the observed Et- As explained in the introduction, the rapidity intervals for particles which 
suppress if> and for those which produce the observed Et do not overlap. Therefore the observables 
mo and Et should not be correlated at all, unless there is cross talk between rapidity intervals. String 
formation and decay is a possible mechanism which leads to correlations between produced particles 
in different rapidity intervals. 

The following calculation of cross talk is based on string formation and their decay within the 
dual parton model 0] for particle production (Fig. 1): Two protons, Pi and P2, with rapidities 
y = and y = Y, respectively, interact via color exchange. After the interaction, two strings form: 
string 1 between the quark with y\ and the diquark of P2 (assumed to have rapidity Y) and string 
2 between the quark with y 2 from P 2 and the diquark from Pi. The probability distribution for the 
occurrence of string i is denoted by Wi(yi), i = 1, 2 and is normalized to 1. The strings fragment into 
hadrons where the number of produced hadrons per rapidity interval is roughly independent of y. In 



5 



a NN collision, a hadron with rapidity y a can arise from each string, provided it covers the rapidity 
y a . Therefore the probability to find a hadron with rapidity y a in the event shown in Fig. 1 is 

P{Va) = j dy x dy 2 Wx{y x ) w 2 (y 2 ) [6(y 2 - y a ) + d{y a - y x )) , (17) 

where each ^-function refers to the contribution of one string. Furthermore, the probability to observe 
two hadrons, one at rapidity y a and another of y\, is 

P(ya,Vb) = Jdy 1 dy 2 w 1 (y 1 )w 2 (y 2 ) (18) 
x [% 2 - Va) 0(y 2 - y b ) + 6(y a - yi) 0(y b - yx)\ ■ 

Let us identify y a with the rapidity of the charmonium y^ and y b with the rapidity yE T of a hadron 
contributing to the observed Et- The probability a x for cross talk (a x < 1) is then given by the 
ratio 

a x {y^ y ET ) = pjjpj > ( 19 ) 

where P(y^,yE T ) gives the probability that a hadron at y^ (comoving with the ip) and a hadron at 
yE T (contributing to Et) are correlated and where P(y^) gives the probability to find a comoving 
hadron at y^ without any further condition. 

The probability distribution Wi(y\) for the occurrence of the string i is taken proportional to the 
quark distribution function fi(xi) via dxifi(xi) = dyiW{(yi), where Xi is the fractional momentum of 
the quark. If we confine ourselves to the low values of x^ then fi(xi) ~ x i 2 and 

wi(yi) = ce~* yi , 

w 2 (y 2 ) = ce-^ Y - y2 \ (20) 

where the normalization drops in the ratio eq. (0). In the experiment under consideration, we are 
interested in the cross talk between the comoving hadrons (3 < y^ < 4) and the hadrons in transverse 
energy E T (1 < y Er < 2.3). Instead of integrating a x over these intervals, we evaluate eq. ( |i~9l) at the 
mean values y^ = 3.5, yE T = 1.65. A straightforward calculation then leads to 

a x (3.5, 1.65) = 0.82. (21) 

The high value of the coefficient for cross talk (a x = 1 corresponds to complete cross talk) reflects 
the property of the probability distributions Wi in that they are largest for those strings which span 
the full rapidity range. In our calculation, we have neglected the strings involving sea quarks. Since 
these strings are shorter, their inclusion would reduce the value of a x . 



With this result, the factor e(E T ) from eq. (|Tj) has to be modified to 

E T {b) 



e(E T ) = l + a .a x ( Et ~ (22) 
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with «o from eq. flI5|) and ot x from eq. (TT9|). Eq. is the final result of our paper. 

We now look at the experimental J/ip suppression as a function of Et for Pb + Pb collisions at 
158 A GeV and compare with it the results of several calculations in order to show the importance 
of the effects discussed in this paper (Fig. 2). We have repeated the calculation of [0 with the form 
given in eq. (^) for the anomalous suppression. However, the expression e = Et/Et in front of n p (b, s) 
is replaced by the expression e derived in eq. (|22| ) with a = «o • ct x which describes the influence of 
the shape of the probability distributions and of the incomplete cross talk. Then a = means that 
no fluctuations are taken into account while a = 1 gives the maximal influence of fluctuations as in 
H [7|. A reasonable estimate may be a = 0.8, but we also give a curve for a = 0.4. While a = 1 and 
a = 0.8 are compatible with the data, a = 0.4 and smaller values are definitely ruled out, provided 
the effect discussed in ref. || is not too important. 

In this paper we have derived and discussed the phenomenological prescriptions by Capella et 
al. H and Blaizot et al. |7j to calculate J/ip suppression in the region of large values of Et where 
the transverse energy fluctuates. The derivation has shown which physical parameters determine 
the strength with which fluctuations in (Ex — Et)/Et influence J/ip suppression. These are the 
width parameters for the multiplicity distributions of comovers and for the Et distributions and the 
correlation function for incomplete cross talk. We conclude that the ad hoc factor, eq. ([!]), holds to 
a good approximation. 
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Figure 1: Particle production for proton-proton collisions within the dual parton model |TT 
protons, Pi and P2 



Two 



with rapidities y = and y = Y collide and exchange color after which two 
strings form. We study production of hadrons at the rapidities (which comove with the ip) and 
yE T (which contribute to the transverse energy E T ). 
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Figure 2: The experimental values for the ratio J/ip/DY in Pb + Pb collisions at 158 A GeV [[[]]. 
The curves are calculated like in the paper by Blaizot et al. [0 and differ by the degree a to which 
fluctuations are taken into account a = «o • a x, eq. (|2l[) . a = : Fluctuations are not accounted 
for. a = 1 : Fluctuations influence suppression fully, a = 0.8 and a = 0.4 : Partial reduction of the 
importance of fluctuations. 
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